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In these lectures I explain how chiral symmetry of continuum QCD naturally leads 
to a class of lattice regularisations known as twisted mass QCD (tmQCD). As 
compared to standard Wilson quarks, its advantages are the absence of unphysi- 
cal zero modes, the possibility to circumvent lattice renormalisation problems and 
automatic 0(a) improvement. On the other hand, the physical parity and flavour 
symmetries are explicitly broken. I discuss these aspects and then turn to the the- 
ory in a finite space-time volume with Schrodinger functional boundary conditions. 
Again, chiral transformations of the continuum theory may be used as a guide to 
formulate an alternative lattice regularisation of the Schrodinger functional, with 
interesting applications to renormalization problems in QCD. 



1. Introduction 

In recent years, twisted mass QCD (tmQCD) has become a popular variant 
of lattice QCD with Wilson-type quarks Initially designed to render 

the (partially) quenched approximation well-defined through the elimina- 
tion of unphysical zero modes, it was soon realised that tmQCD could also 
be used to circumvent some notorious lattice renormalization problems 
Later, Frezzotti and Rossi observed that scaling violations in tmQCD can 
be reduced to O(a^) without the need for all the 0(a) counterterms required 
with standard Wilson quarks (a being the lattice spacing). This property, 
referred to as "automatic 0(a) improvement", has attracted further atten- 
tion and a number of groups have started large scale numerical simulations 
using tmQCD. In these lectures I do not attempt to review this work in 
progress''. Here I would rather like to give an introduction to the basic con- 



* based on lectures given at the school "Perspectives in Lattice Gauge Theories", Nara, 
Japan, October 31- November 14, 2005 
''See ^ for a review and further references. 
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cepts. This includes in particular a discussion of 0(a) improvement and 
the question whether it is compromised by currently used non-perturbative 
renormalization procedures based on the QCD Schrodinger functional (SF- 
schemes). In fact, the standard Schrodinger functional boundary conditions 
turn out to be difficult to reconcile with automatic 0(a) improvement and 
the construction of an alternative set-up for the Schrodinger functional may 
therefore be advantageous. 

This writeup is organised as follows: I start with the interplay between 
the choice of the quark mass term and the form taken by parity, flavour 
and chiral symmetry transformations (sect. 2). After a reminder of stan- 
dard Wilson quarks and the problem of unphysical zero modes (section 3), 
lattice tmQCD is introduced in sect. 4. Based on the formal continuum the- 
ory a dictionary between tmQCD and QCD correlation functions is readily 
established, which is expected to hold between properly renormalised cor- 
relation functions. It then becomes clear how to by-pass certain renormal- 
ization problems of standard Wilson quarks (sect. 5), and the computation 
of Bk is discussed in some detail. In sect. 6 automatic 0(a) improvement 
of tmQCD is analysed using Symanzik's effective theory. Potential prob- 
lems of tmQCD associated with flavour and parity breaking arc shortly 
mentioned in sect. 7. In sect. 8, the properties of Schrodinger functional 
renormalisation schemes (SF schemes) are discussed. Motivated by the 
clash of the standard set-up with automatic 0(a) improvement and by the 
slow dccoiipling of heavy quarks in mass-dependent SF schemes, a modi- 
fied definition of the Schrodinger functional is proposed, and its effective- 
ness regarding 0(a) improvement is illustrated in an example taken from 
perturbation theory. Section 9 contains some conclusions. 

2. Continuum QCD and chiral transformations 

Let us consider the continuum action of QCD with N{ = 2 massless quarks''. 
Decomposing the action into a pure gauge and a fermionic part, S = Sg+Sf, 
we here focus on the fermionic part. 



The quark and antiquark fields xp,ip are flavour doublets, interacting min- 
imally with the gluon field via the covariant derivative = dfj, + A^^. 



''Conventions used for Euclidean 7-matrices in 4 dimensions: {7^,71/} = 2Spti,, 7J = Ji^i, 
where /i, 1/ = 0, 1,2,3, and 75 = 707l7273, cr^i^ = fb^iTL-]- 





(1) 
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The massless fermionic action has a global chiral-flavour SU(2)xSU(2) in- 
variance, corresponding to the transformations, 

V' ^ V' = exp(iw^T"/2) exp(iw^75TV2)i'. 

= ,/;cxp(7;cj^75tV2) cxp(-it^^^T'72), (2) 

where r" (a = 1,2,3) are Pauli matrices, are transformation param- 
eters'^. This notation distinguishes the axial from the vector generators 
(corresponding to the flavour or isospin SU(2) subgroup) in a standard 
way. 

A quark mass term breaks the chiral flavour symmetry explicitly, leaving 
only the vector or isospin symmetry intact. The above notation for the 
symmetry transformations was introduced with the standard quark mass 
term tl)^ in mind, but e.g. the choice 

= cxp(«Ci;^75t")?/; = cos (w^t )■(/>?/' + is\n{LL)jC)u'\'tp^ir,T°'i). (3) 

would be completely equivalent. Here, uja denotes the modulus of 
{uj\,io\,ijO^ and u\ = u:\lu: a is a unit vector. In fact, it is only after 
the introduction of the quark mass term that the distinction between ax- 
ial and vector symmetries acquires a meaning. By definition, the vector 
symmetry transformations are those which leave the quark mass term in- 
variant. Similarly, the quark mass term is supposed to be invariant under 
parity transformations. As a consequence, the form of a symmetry trans- 
formation depends on the choice of the mass term. While a standard mass 
term implies that a parity transformation can be realised as 

■'/'(a;o,x) ^ 7o^!>(xo, -x), ^/^(ajojx) ^ !/)(2;o, -x)7o, (4) 

the alternative choice of (3) for the mass term means that a parity trans- 
formation will look more complicated, for instance 

V'(a;o,x) 7oexp(iw^75r")V'(a;o, -x), 

Vi(a;o, x) 'ijj{xo, -x) exp(iw^75r")7o. (5) 

Similarly, the isospin transformation obtained with a standard mass term 
corresponds to (2) with all axial transformation parameters set to zero, 
coa = (whence the notation), whereas the mass term (3) leads to the 
much less intuitive formula 

ij exp(-iw^75r"/2)exp(ia;^rV2)exp(ia;^75r72)V', 

Vi ^ Vi exp(ia;^75T"/2) e^p{-iuj^T'' /2) exp(-ja;^75r72), (6) 



'^Summation over repeated indices a,b = 1, 2, 3 is understood. 
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where Uy {b — 1,2,3) are transformation parameters while uj'\ (a = 1, 2, 3) 
are again fixed. The situation is reminiscent of the choice of a coordinate 
system, and our intuition about the form of symmetry transformations is 
thus based on a particular choice of "field coordinates". Of course, this 
raises the question why one should deviate from the standard choice of 
the mass term. In the continuum and for regularisations preserving chiral 
symmetry there is indeed no point in introducing a twisted mass term, 
for any non-standard choice could be brought into the standard form by 
using an axial rotation, which, being a symmetry of the masslcss theory, 
has no further effects. The situation is different in regularisations which 
break chiral symmetry, such as lattice regularisations with Wilson type 
quarks. One may then obtain different regularisations of QCD which have 
equivalent continuum limits but differ at the cutoff level. This will be made 
more precise a bit later. 

3. Standard Wilson qucirks 

Standard Wilson quarks are characterised by the fermionic lattice action, 

Sf = a''J2i'i^)iDw+mo)^{x), (7) 

X 

3 

^w = E{5(v^ + v;)7M-av;v^}. (8) 

Here, mo is a bare mass parameter and the covariant lattice derivatives in 
the Wilson-Dirac operator are defined as usual (see ref. ^ for unexplained no- 
tation). Assuming A^f quark flavours the lattice action has an exact U(Af) 
vector symmetry, and is invariant under axis permutations, reflections such 
as parity and charge conjugation. Furthermore, unitarity of lattice QCD 
with Wilson quarks has been rigorously established These nice proper- 
ties of standard Wilson quarks come with a price: all axial symmetries are 
explicitly broken by the last term in eq. (8), called the Wilson term. This 
has a number of consequences: 

(1) Linear mass divergence: the quark mass term is not protected 
against additive renormalization, i.e. any renormalized quark mass 
is of the form mR ~ Zminio — Wcr), where the critical mass is lin- 
early divergent, i.e. rucr oc 1/a. 

(2) Axial current renormalization: since axial transformations are not 
an exact symmetry, there is no exact current algebra, and the non- 
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singlet axial current requires a non-trivial multiplicative renormal- 
ization to restore current algebra up to 0(a) effects. 

(3) Definition of the chiral condensate as expectation value of a local 
operator: the renormalised iso-singlet scalar density has the struc- 
ture, 

{tPi>)^ = Zso{^^ + csa-^}. (9) 

In a regularisation which respects chiral symmetry, the additive 
renormalization constant cs would be proportional to am, with 
m being a multiplicativcly rcnormalisable bare quark mass. This 
means that the chiral condensate is well-defined in the chiral limit 
once its multiplicative renormalisation has been carried out. In con- 
trast, with Wilson quarks one first needs to subtract the cubic power 
divergence, even in the chiral limit. 

(4) Cutoff effects: the leading cutoff effects with Wilson-type fermions 
are proportional to a, rather than a^. Again, this is a consequence 
of chiral symmetry breaking. This is easily seen by looking at the 
structure of the counterterms which are to be included for the on- 
shell 0(a) improvement of the theory a la Symanzik ^. 

From a field theoretical point of view this illustrates the proliferation of 
additional counterterms in a case where the regularisation breaks a con- 
timuun symmetry. One should note, however, that there is no remaining 
theoretical or conceptual problem. 

3.1. Wilson quarks and unphysical fermionic zero modes 

Nevertheless, technical problems may arise within the current practice of 
numerical simulations with Wilson- type quarks. This is related to the fact 
that, for a given gauge background field, the massive Wilson-Dirac oper- 
ator Dw + 'TiQ is not protected against zero modes unless the bare mass 
parameter toq is positive. However, due to additive quark mass renor- 
malisation, the masses of the light quarks typically correspond to negative 
bare mass parameters, which leaves the Wilson-Dirac operator unprotected 
against zero modes in the physically interesting region. These modes are 
considered unphysical, since one expects from the continuum theory that 
any non-zero value of the renormalised quark mass prohibits zero modes of 
the Dirac operator. 

It is instructive to look at a typical fermionic correlation function, such 
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as the pion propagator given by 

G-\x -y) = - (i^{x)i/-^i>{x)i,{y)^J-^i^{y)^ 

= J Z)[[/,^,^]e-^^(a;)75YV(a;)V^(y)75YV'(y),(10) 

where r", a = 1,2,3, are the Pauli matrices acting in flavour space and 
Z = (1). It is convenient to introduce the operator, 

Q = 75(r'w + mo), Q = Q\ (11) 

which acts in single flavour space. Integrating over the quark and anti-quark 
fields one obtains 

G^'ix -y)= ^-S'^'Z-^ j D[U]e-^<' det (Q^) tr [Q-\x,y)Q-\x,y)\ , 

(12) 

where the flavour structure has been reduced analytically and the remaining 
trace is over colour and spin indices. The important point to notice is that 
the resulting expression is never singular. Denoting the eigenfunction of Q 
for a given eigenvalue A, by (pi{x), the pion propagator takes the form 

X E ^7'K'^ji^)v*jiy)'pk{xMiy)- (i3) 

In other words, the eigenvalues in the denominator are always compensated 
by corresponding factors from the determinant. The limit of vanishing 
eigenvalues is always regular and a strict lower bound on the eigenvalue 
spectrum is not required for the theory to be well-defined. 

However, the absence of a lower bound on |Aj| may still lead to technical 
problems, either due to the use of unphysical approximations or due to the 
set-up of numerical simulations: 



3.1.1. Quenched and partially quenehed approximations 

As the computational cost for the generation of a an ensemble of gauge 
field configurations is dominated by the inclusion of the quark determi- 
nant, a widely used approximation consists in omitting the determinant 
when taking the average over gauge fields. The quark propagators with the 
eigenvalues in the denominator may then become singular, and gauge field 
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configurations wtiere this happens are called "exceptional". The example 
in figure 1 taken from * shows the ensemble average of the pion propagator 
over all gauge configurations but a a single exceptional one (dashed line), 
where the propagator deviates dramatically from the average (dots and 
solid line). The inclusion of the exceptional configuration in the ensemble 
average would lead to much larger errors, while its omission invalidates the 
Monte Carlo procedure. In principle one should say that the quenched ap- 
proximation with Wilson type quarks is ill-defined, since zero modes are 
bound to occur if the ensemble of gauge configurations is large enough. 
However, the frequency of near zero modes depends very sensitively upon 
the bare quark mass and is in fact a function of the lattice size and all the 
other bare parameters in the lattice action. One may therefore think of the 
quenched approximation as being operationally defined, if for an ensemble 
of, say, a few hundred configurations the problem is typically absent. "Safe" 
parameter ranges may then be quoted for a given action, but this situation 
is clearly imsatisfactory. In particular, as the problem is not sharply de- 
fined, one may always be unlucky and encounter near zero modes even at 
parameter values which have previously been considered safe. In practice 
it is this problem which has limited the approach to the chiral limit, rather 
than finite volume effects due to the pions becoming too light. 

Obviously, the problem is expected to disappear once the quark determi- 
nant is properly included. Usually this is done by including the complete 
determinant in the effective gauge field measure used in the importance 
sampling, and the probability for a gauge configuration to be included in 
the ensemble becomes proportional to the eigenvalues. Exceptional con- 
figurations are then never produced. However, even in this case, one is 
often interested in varying the valence quark masses independently of the 
sea quark masses, a situation which is referred to as the partially quenched 
approximation. One may also have different numbers of valence and sea 
quarks, or Wilson valence quarks and sea quarks of a different kind. In all 
these cases one expects similar problems with unphysical zero modes as in 
the quenched approximation. 

3.1.2. Potential problems in the Hybrid Monte Carlo algorithm 

Most numerical simulations use some variant of the Hybrid Monte Carlo 
algorithm ^. Integrating the molecular dynamics trajectories in fictitious 
phase space then requires the evaluation of the fcrmionic force term and 
thus the inversion of the Dirac operator at each step in molecular dynamics 



February 7, 2008 10:0 Proceedings Trim Size: 9in x 6in 



ms 



8 




10 20 30 40 50 60 



t 

Figure 1. The pion propagator vs. time separation from a quenched simulation on a 
32^ X 64 lattice at ^8 = 6.2 (cf. text for an explanation). 

time. The force term may become very large if an exceptional configuration 
is encountered, and the molecular dynamics integrator tends to become un- 
stable if the product of the force and the step size exceeds a certain critical 
value To avoid this situation one may hence be forced to decrease the 
step size to very small values thereby increasing the cost of the simulation. 
It is likely that this problem was at the heart of the difficulties encountered 
in the past with simulations of Wilson type quarks However, various de- 
velopments over the past few years seem to have solved this problem (see 
for a recent account of current simulation algorithms and cost estimates and 
for further discussion). 

4. Twisted mass lattice QCD 

Initially the main motivation for introducing a twisted mass term was the 
problem with zero modes discussed in the previous section. The lattice 
action for a doublet ip oi N{ = 2 mass degenerate quarks is now given by 



Sf = a^y^V>(a;)(J>vy -|- mo -|- i/Uq75T^)V'(a;), 

X 



(14) 
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where /iq denotes the bare twisted mass parameter. It is easy to see that 
the presence of this parameter eUmiates any unphysical zero modes, for 

det (Dw + mo + iMq75''"^) = det 

= det(Q2 + > 0. (15) 

The difference in the determinant already shows that twisted mass and 
standard QCD cannot be the same regularisation. In fact, any attempt 
to perform an axial rotation so as to eliminate the twisted mass term will 
rotate the Wilson term in eq. (8), too. the equivalence between both regu- 
larisations can therefore only be expected to hold in the continuum limit. 
We will discuss this more in detail below. 

Here it suffices to say that the chiral flavour symmetry of twisted mass 
QCD is reduced to an exact U(l) symmetry with generator Further- 
more, charge conjugation, axis permutations and reflections combined with 
a flavour permutation, e.g. 

■^(a;o,x) ^7orV(a;o,-x), V;(a;o,x) -> i^(a;o, -x)7ot\ (16) 

arc exact symmetries. Finally, the construction of a positive and self-adjoint 
transfer matrix for standard Wilson quarks can be generalised to twisted 
mass QCD, provided /Zq is real and the usual condition on the standard 
bare mass parameter, |«;| < 1/6, with k = (2amo + 8)~^ is satisfied ^. 

4.1. Equivalence between tmQCD and QCD 

Taking the continuum limit, we see that the fermionic continuum action of 
tmQCD, 

Sf = J dS ^Pix){I^ +m + i/iq75T3)V'(a^), (17) 

can be related to the standard action by a global chiral field rotation, 

V'' = i?(a)V', ^l)' = ^R{a), i?(a) =exp (^?:a75y^ . (18) 

Choosing the angle a such that tana = jJ-q/m, the action for the primed 
fields takes the standard form, 

S'f = J d*x ^j'{x){ip+M)^\x), M = ^m? + nl- (19) 

In QCD all physical obscrvablcs can be extracted from gauge invariant cor- 
relation functions of composite fields. We would therefore like to study the 
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relationship between correlation functions in tmQCD and standard QCD. 
To this end we introduce polar mass coordinates, 

m = Mcos(Q:), fiq = Ms'm{a), (20) 

and consider the correlation functions labelled by (M, a), 

(0[V'.V^])(M,a) = ^"'/ D[U,i^,i^] 0[^,^] e-^I-'^-l. (21) 

Treating the functional integral like an ordinary integral we change the 

variables to tjj' and -ip' of eq. (18) and re-labcl these new integration variables 
to ^ and tjj afterwards. In this way wc arrive at the identity, 

Vi])(M,o) = mc^)])^M,c) ■ (22) 

To go a step further, we now assume that the functional 0[tjj, ijj] consists of 
factors which are members of a chiral multiplet. Considering such a field 
in the representation r, the transformation of ip and -ip by R{a) 
induces the transformation of by R^^^ (a) in the representation r, 

[i?(a)^, i,R{a)] = 41(a)4^^ [V, V^]- (23) 
For n-point functions of such fields, one obtains the identity, 



^M,0) = n<k(") (<'^---'^B:^)(M.a)- (24) 



\i=l 



The correlation fnnctions in standard QCD labelled by (M, 0) are just linear 
combinations of those in twisted mass QCD, labelled by (M, a). The inverse 
relation can be obtained by inverting the matrices i?('''(a). This is trivial, 
as the axial rotation (18) forms an abelian subgroup of the chiral flavour 
group, so that [R^'''\a)]^'^ = R^'^^—a). Examples of chiral multiplcts arc 
the non-singlet currents (A^, V"^) or the non-singlet axial density combined 
with the singlet scalar density, {^S^,P°'). In terms of quark fields one has 

and one may then easily infer the transformation behaviour of these chiral 
multiplets: 

A'^ = cAl + sV^, = cV^ + sAl , 

A'^=cAl-sV^, V;,' = cV^ - sAl, 

K=Al, V;,' = V^, (26) 

pm^pa^ (a =1,2), S'° = cS° + 2isP^, 

p/3 ^ ^p3 + ^glsO_ 
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Here the notation O' = 0\^' ,'^'\, c = cos(a), s = sin(Q!) was used. For a 
correlator of A^^{x) and P^{y) in standard QCD this means 

<4(x)pi(y))(^_„^ =cos(«) {Al{x)P\y))^^^^^ 

+ sm{a) {V^{x)P\y))^^^^^. (27) 

In other words, eqs. (26) relate an insertion of the primed fields into stan- 
dard QCD correlators to the insertion of the corresponding r.h.s. into 
tmQCD correlators. In particular, we note that the PCAC and PCVC 
relations in the physical basis 

d^A'- = 2MP'% = o, (28) 

are equivalent to linear combinations of their twisted counterparts, 
d^Al = 2mP» + S^'^ifigS^, 

d^VH = -2/Xge3»''P\ (29) 

In conclusion, the formal continuum theory provides us with a dictionary 

between correlation functions in standard and twisted mass QCD. However, 
all these considerations have been quite formal, and we need to specify how 
such a dictionary carries over to the renormalized theories. 

4.2. Beyond the formal continuum theory 

To clarify this question let us suppose that tmQCD is regularised on the 
lattice with Ginsparg- Wilson quarks, where chiral and flavour symmetries 
are the same as in the continuum. Identities such as Eqs. (24) may then 
be derived in the hare theory. If in addition, we start from a finite vol- 
ume with, say, periodic boundary conditions for all fields, the functional 
integral becomes a finite dimensional Grassmann integral. Therefore, these 
identities are no longer formal, but on firm mathematical grounds, and 
all one has to show is that the renormalisation procedure can be carried 
out such that they continue to hold in the renormaliscd theory. This is 
straightforward, as one just has to make sure that all members of a given 
chiral multiplet are renormaliscd in the same way, and that the multiplica- 
tive renormalization constants do not depend on the twist angle a. This 
can be achieved e.g. by imposing renormalisation conditions in the massless 
limit. Hence, in this case, the dictionary introduced above holds between 
the renormaliscd correlation functions of both theories. Assuming univer- 
sality to hold beyond perturbation theory, this establishes the equivalence 
of both versions of QCD at the non-perturbative level, since any other 
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regularisation, chirally symmetric or not, will then lead to the same renor- 
malised correlation functions up to cutoff effects. While there is no reason 
to doubt that universality holds generally, one should be aware that it has 
rigorously been established only in perturbation theory and for selected 
regularisations (e.g. lattice regularisations with Wilson type quarks ^^). 

4.3. Lattice tmQCD with Wilson quarks 

In tmQCD on the lattice with Wilson quarks the axial transformation re- 
lating continuum tmQCD to standard QCD is not an exact symmetry. 
Therefore, equivalence can only be expected to hold in the continuum limit, 
i.e. for properly rcnormalized correlation functions and up to cutoff effects. 
The lattice symmetries imply the counterterm structure, with the following 
result for the renormalised parameters, 



It is a priori not obvious how the twist angle a should be defined from 
the mass parameters. The key observation is that chiral symmetry can 
be restored in the bare lattice theory up to cutoff effects, by imposing ax- 
ial Ward identities as normalisation conditions . This fixes the relative 
renormalization of all members of a chiral multiplct, such as Z^/Zy for the 
symmetry currents'^, or Zgo/Zp for the iso-triplct axial and the iso-inglet 
scalar densities. Note that such ratios arc scale independent functions of qq 
only, which are expected to converge to 1 in the continuum limit with a rate 
(7q oc — 1/lna. In particular, these ratios do not depend upon the quark 
mass parameters and may therefore be determined in the massless limit 
where the tmQCD and standard QCD actions coincide. The connection 
between the mass parameters and chiral Ward identities is established by 
choosing renormalisation schemes such that the PCAC and PCVC rela- 
tions hold, with the renormalised currents and and axial density, and the 
renormalised mass parameters. The renormalization constants may then 
be shown to satisfy the identities, Z^ = Z^q and Z^ = Z^^. With these 
conventions it is clear that the ratio of renormalised mass parameters is 
known once the critical mass and the ratio Zgo /Zp are given. 



•J. ry '2, 

9r — ^sffo ' 



rriR = Zm{mo - rricr) 



(30) 



tana = = 



Zp mo — m, 



'cr 



(31) 



'^Zy = 1 only holds if the (partially) conserved point-split vector current is used. 
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Besides the ratio of renormalization eonstants one thus needs to determine 
the critical mass. In practice this can be done by measuring a bare PCAC 
mass m from correlation functions with some external field O, 

(dr,A](x)0) 

•"'-wikw' '''' 

and by using the relation 

rriR = Zp^Zpja, m = Z^^Z„Zp(mo - rricr). (33) 

Alternatively one may use the measured bare PCAC quark mass m to 
obtain a directly, 

tana = /iq/(ZAm), (34) 

provided one has previously determined Za- Already at this point one 
notes that the choice a = 7r/2 is special, as in this case one merely needs 
to determine the critical mass. The choice a = 7r/2 is referred to a full or 
maximal twist, because the physcial quark mass is then entirely defined by 
the twisted mass parameter /Xq. 

Having determined the twist angle, and the relative renormalizations 
within chiral multiplets, chiral symmetry is restored up to cutofl[ effects for 
the correlation functions of members of these multiplets. In a second step 
one just needs to make sure that this property of the bare theory is not 
compromised by the renormalization procedure, i.e. one is in a similar situ- 
ation as in the bare theory with Ginsparg- Wilson quarks. Proceeding in the 
same way, the formal identities of subsect. 4.1 will hold in the renormalised 
theory. 

An important point to notice is that the twist angle a is a new parameter 

which reflects the freedom to choose a direction in chiral flavour space for 
the explicit chiral flavour symmetry breaking. Our physical interpretation 
is such that by definition only the axial generators are broken by the mass 
term thus defining the residual vector symmetry. With Wilson quarks at 
non-zero a there is an additional breaking of flavour symmetry by the 
Wilson term, which is expected to disappear in the continuum limit, just 
like chiral symmetry is restored with standard Wilson quarks. In order to 
define the continuum limit properly one must make sure that cutoff effects 
are a smooth function of /3 = 6/50 ■ In general this can be achieved by taking 
the continuum limit at constant physical conditions. For instance one may 
keep TOtt/J^tt constant as /? is varied. However, in tmQCD this observable is 
a function of two mass parameters, or, equivalently of one mass parameter 
and the twist angle. It is crucial that the twist angle is kept constant as 
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the continuum limit is taken, since the twist angle labels different lattice 
regularisations of two-flavour QCD. In particular, if a is changed from 
one /3- value to the next, there is no reason to expect a smooth continuum 
approach and a continuum extrapolation may become impossible. 

5. A few applications of tmQCD 

The relations between tmQCD and standard QCD correlation functions can 

be used to by-pass certain lattice renormalization problems of standard Wil- 
son quarks (cf. sect. 3). As the different operators of a continuum chiral 
multiplet are not necessarily related by lattice symmetries, their renor- 
malisation properties can be very different. Moreover, the renormalisation 
properties do not change in the presence of (twisted or non-twisted) mass 
terms except when power divergences are present. Excluding these cases it 
is thus sufficient to renormalise a given composite field in the chiral limit 
where the actions of tmQCD and standard Wilson quarks coincide. One 
may then choose the operator with the best renormalisation properties that 
can be related to the desired standard QCD operator by the dictionary es- 
tablished earlier. Moreover, it may not even be necessary to match the 
operators directly. In principle, it is enough to match the desired correla- 
tion function up to cutoff effects. Perhaps these remarks become clearer by 
going through a few examples: 

5.1. Computation of F^^ 

Both the pion mass m,r and the pion decay constant can be obtained 
from the long distance behaviour of the 2-point function 

((AR)J(x)(PR)i(2/)a)(^^_g) =cos(a) {{A^)l{x){P^)Hy)) ^^^^^^ 

+ sin(a) (vS{x){P^)\y)) .(35) 

\ / (Mr,c<) 

The problem with the standard Wilson computation on the l.h.s. is that 

the axial current requires a non-trivial renormalisation, which needs to be 
determined from Ward identities, as done e.g. in On the other hand the 
vector current is protected against such a rescaling since it is conserved 
at /Uq = 0. At a = 7r/2 the axial current is mapped to the vector current and 
one may thus avoid the current renormalisation by computing the vector 
correlation function in tmQCD. It is in fact not necessary to set a = 7r/2; 



February 7, 2008 10:0 Proceedings Trim Size: 9in x 6in 



ms 



15 

when inverting the relation (35), 

{{Vo\^){PK)\y)),^^ =cos(a) {{Vo\^)iPK)\y)),^^ 

\ / (Mh,q) \ / (Mr,0) 

+ sin(a) ((AR)i(x)(PR)i(y))(^^Q). (36) 

one notices that the first term on the r.h.s violates both parity and flavour 
symmetry of standard QCD. On the lattice this correlation function there- 
fore contributes at most an 0(a) effect. One may thus obtain i<V at values 
a ^ tt/2 by computing the l.h.s of this equation. Finally, it should be 
mentioned that the exact PCVC lattice relation, 

d;v; = -2Mqe'«^P^ (37) 

may be used to replace the vector current by the axial density. Summing 

over X, translation invariance eliminates the spatial part of the divergence, 
and the time derivative reduces to a multiplication by TOtt at large time 
separations The results of a quenched computation along these lines 
are shown in figures 2 and 5. 

5.2. Direct determination of the chiral condensate 

A computation of the chiral condensate from the local scalar density has 
never been performed with Wilson quarks, due to the cubic divergence (9) 
which persists in the chiral limit. In tmQCD the role of the scalar density 
is played by the axial density, i.e. one expects the relation 

(38) 

Again, the first term on the r.h.s. vanishes up to 0(a) due to parity, so that 
the computation of the l.h.s. yields the chiral condensate up to the factor 
(— j/2) sin(a). This is advantageous as the renormalised axial density is of 
the form, 

(Pnf = Zp (P3 + /,qcp a-^) , (39) 

i.e. the power divergence vanishes for /itq = 0. Still, in order to determine 
the condensate one needs to perform first the infinite volume limit followed 
by the fj,q = 0, mg = m„ limits (at fixed a) and the continuum limit, which 
remains a rather delicate task. In particular, the chiral limit is complicated 
by the fact that the uncertainty in sin(a) increases as the quark mass is 
decreased, due to the intrinsic 0(a) ambiguity of rrici. In practice this 
means that one has to extrapolate to the chiral limit from some distance, 
but this is anyway required for finite volume effects to remain small. 
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Figure 2. Quenched continuum results for F-n. The plot also illustrates the absence of 
the zero mode problem in tmQCD, as much smaller pion masses could be reached than 
with standard 0{a) improved Wilson quarks. 



5.3. The computation of Bk 

Four-quark operators provide an interesting playground for mappings be- 
tween tmQCD and standard QCD. We start with the Bk parameter which 
is defined in QCD with dynamical m, d, s quarks by, 



<^(V-A)(V-A) 



The local operator 



'^(V-A)(V-A) 



(40) 



(41) 



is the effective local interaction induced by integrating out the massive 
gauge bosons and t-^h- and c-quarks in the Standard Model. The transi- 
tion between the pseudoscalar states and does not change parity. 
Therefore, only the parity-even part in the effective operator. 



0(V-A)(V-A) = OvV+AA — OvA+AV , 

parity-even parity-odd 



(42) 
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contributes to Bk- With Wilson type quarks, the operators Oyv+AA and 
OvA+Av are renormahsed as follows 

4 

(Ovv+aa)r = ^vv+aa|Ovv+aa + ^2 ^» (^^) 

(Ova+av)r — ^va+avOva+av- (44) 

While the parity-even component mixes with four other operators of di- 
mension 6, the parity-odd component only requires multiplicative renor- 
malisation, due to CP and flavour exchange symmetries This raises the 
question if one can by-pass the mixing problem by exchanging the roles 
of both operators through the introduction of twisted mass terms. This is 
indeed possible, but one first needs to introduce the strange quark. The 
simplest possibility consists in adding a standard s-q\iark to a twisted quark 
doublet 7p of the light up and down quarks, which are thus taken to be de- 
generate. The corresponding continuum Lagrangian is given by 

/: = -ip{Ip+m + iiig-fr,T^)ilj + sil^+m,)s, (45) 

and, passing to the physical basis of primed fields, one finds 

Oyv+AA = cos(a)Ovv+AA - * sin(a)OvA+Av 

= -iOvA+Av {a = 7r/2). (46) 

At full twist, we thus get a direct mapping between both operators, 
i.e. OvA+Av in twisted mass QCD at 7r/2 is interpreted as Ovv+aa in 
standard QCD. A second possibility consists in exchanging the roles of up 
and strange quark, i.e. one considers a twisted doublet of strange and down 
quarks and a standard u-quark. In this case one finds 

O'w+AA = cos(2a)Ovv+AA - i sin(2a)0vA+Av 

= -iOvA+Av {a = 7r/4), (47) 

i.e. the same mapping is obtained, but with the twist angle a = 7r/4. Sev- 
eral comments are in order: while both options, referred to as 7r/2 and 
7r/4 scenarios respectively, are possible, the second one is clearly more re- 
mote from reality, as it assumes mass degenerate down and strange quarks. 
However, this is precisely the limit in which most lattice calciilations to 
date have been performed. The justification rests on chiral perturbation 
theory where a a weak dependence upon the strange-down mass difference 
is predicted. Moreover, in the quenched approximation, any deviation from 
the degenerate case leads to an unphysical logarithmic quark mass depen- 
dence 



February 7, 2008 10:0 Proceedings Trim Size: 9in x 6in 



ms 



18 

5.3.1. Renormalisation of OyA+AV 

Whatever the chosen strategy, the operator which requires renormalisation 
is OvA+AV- The renormalisation is multiplicative, and the general strat- 
egy of ^'^ can be applied. The scale evolution of the operator in a few 
Schrodinger functional schemes has been traced in the quenched approxi- 
mation over a wide range of scales (for first results with N[ — 2 sea quarks 
cf. ^''). The result is shown in figure 3. It thus remains to calculate the 
bare matrix element for Bk at various values of /?, and, after multiplica- 
tion with the Z-factor at the low energy scale, perform the continuum limit 
extrapolation. In the continuum limit one may then use the known scale 
evolution to reach the truly perturbative regime where contact is made with 
the perturbative renormalisation schemes of the continuum. 



SF; scheme, Nf=0 _ 
1-loop 7, 2-loop jS 




Figure 3. The data points show the non-perturbatively computed scale evolution of Bk 
in the SF scheme. Also shown are two perturbative approximations. 
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5.3.2. Results for Bk in the quenched approximation 

Both scenarios have been implemented in the quenched approximation with 
lattice spacings a = 0.05 — 0.1 fm and lattice sizes up to L/a = 32 If 
one sticks to mass degenerate down and strange quarks the 7r/2 scenario re- 
quires some chiral extrapolation, due to the problem with unphysical zero 
modes (recall that the s-quark remains untwisted). In the 7r/4 scenario 
the zero mode problem is eliminated and the kaon mass can be reached 
by interpolation, provided the finite volume effects are small enough. This 
is the case with all lattice spacings except the finest one, where some ex- 
trapolation is required. A combined continuum extrapolation to both data 
sets, linear in a, leaving out the data at the coarsest lattice spacing led to 
the result Bk = 0.789(46) where Bk denotes the renormalisation group 
invariant S-paramctcr. Unfortunately, the twist angle at /? = 6.1 had not 
been tuned precisely enough, a fact that was only noticed after publication 
of A new analysis indicates that higher than linear lattice artefacts are 
still significant at /3 = 6.1. As the data set is not sufficient to fit to both a 
and terms, it was decided to discard the data at /3 = 6.1, too, with the 
result 2^ 

5^ = 0.735(71) o Sf^(2GeV) = 0.534(52), (48) 

which is compatible with the earlier result, albeit with a larger uncertainty. 
In conclusion, the quenched result for Bk has a total error of almost 10 
percent, which includes all systematic effects (renormalisation, chiral inter- 
or extrapolations, continuum extrapolation) except quenching and the fact 
that the valence quarks are mass degenerate. However a variation of the 
mass difference up to (Mg — M^/{Ms + Mj) « 0.5 did not show sizeable 
effects. While the error could still be improved by including data at a 
finer lattice spacing it seems fair to say that further progress requires the 
inclusion of sea quark effects. 

5.4. Further applications 

Twisted mass QCD does not provide a general recipe for by-passing the lat- 
tice specific renormalisation problems of Wilson quarks. Rather, one needs 
to discuss on a case by case basis whether it can be advantageous to use 
some variant of tmQCD. For further applications to four-quark operators 
and K n transitions I refer the reader to 26,27 \Yhile the first reference 
insists on an equal treatment of sea and valence quarks, the second paper 
explores a mixed action approach, where the valence quarks are chirally 
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Figure 4. Quenched lattice data for both scenarios, a = n/2 and a = tt/4. Also shown 
is the combined continuum extrapolation, leaving out the data at the two coarsest lattice 
spacings. 

twisted individually, independently of the sea quark action. This yields a 
much greater flexibility and allows for a complete elimination of lattice spe- 
cific mixings and subtractions, even including 0(a) improvement. Finally, 
similar considerations apply to QCD with static 6-quarks, where the mix- 
ing of four-quark operators is considerably simplified by twisting the light 
quarks (see for a recent review and further references). 

6. 0{a) improvement and tmQCD 

Given that the quenched approximation is currently being overcome, and 
the zero mode problem for algorithms can be alleviated, there remain es- 
sentially two arguments in favour of tmQCD as opposed to standard or 
0(a) improved Wilson quarks: the first consists in the possibility to by- 
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pass renormalisation problems, as explained in the preceding section. The 
second, is the property of "automatic 0(a) improvement" at maximal twist 
(i.e. a = 7r/2), as first observed by Frezzotti and Rossi ^. I will explain 
this point more in detail below, after a brief reminder of the situation with 
standard Wilson quarks. 

6.1. 0(a) improvement of Wilson quarks 

In lattice QCD with Wilson quarks, results are typically affected by 0(a) 
lattice effects, which is to be contrasted with staggered or Ginsparg- Wilson 
quarks where the leading cutoff effects are quadratic in a. As illustrated 
by the Bk determination described above, linear lattice artefacts render 
continuum extrapolations more difficult, and it would be nice to get rid of 
them altogether. This is possible by introducing 0(a) counterterms to the 
action and the composite operators such that 0(a) effects are cancelled in 
on-shell quantities. The basic idea goes back to Symanzik while the 
restriction to on-shell quantities in gauge theories has been first advocated 
by Liischer and Weisz When applied to Wilson quarks it turns out 
that 0(a) improvement of the spectrum (particle masses and energies) can 
be achieved by adding a single counterterm to the action, the so-called 
Sheikholeslami-Wohlert (SW) or clover term, i'il)aij,vF^i,tpj where F^i, is the 
gluon field tensor This term is of dimension five and therefore comes 
with an explicit factor a when included in the lattice action density. 

While 0(a) improvement of the spectral quantities is quite economical, 
one is often interested in matrix elements of composite operators, and each 
operator comes with its own set of 0(a) counterterms, all of which have to 
be tuned in order to cancel the linear lattice artefacts. While this may still 
be possible for quark bilinear operators, the counterterms quickly proliferate 
in the case of 4-quark operators, and 0(a) improvement becomes completely 
impractical if the quarks are taken to be mass non-degenerate (cf. '^^). 

6.2. Automatic 0(a) improvement of tmQCD in a finite 
volume 

The Symanzik effective theory can also be applied to tmQCD and a list of 
0(a) counterterms for the action and a few quark bilinear operators can be 
found in The observation in ^ is, that at maximal twist, all the 0(a) coun- 



"^On the lattice the field tensor is usually discretised using four plaquette terms in the 
(/i, !^)-plane whence the name "clover term" . 
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tertcrms become irrelevant in tlic sense that they can at most contribute at 
O(a^). The argument for automatic 0(a) improvement can be made such 
that it only relies on Symanzik's effective continuum theory 33,5,34,35 
simplify the discussion, let us first assume that the space-time volume is 
finite, so that spontaneous symmetry breaking is excluded and all observ- 
ables are analytic in the quark mass parameters. We furthermore assume 
that we have tuned some PCAC current quark mass mpcAC = 0, i.e. the 
renormalized standard mass parameter vanishes up to 0{a) effects. Then 
Symanzik's effective continuum action is given by 

Ses = So + aSi + 0{a^), So = j d'^x i){ip + iiiq-inT^)^), (49) 

where Sq is the maximally twisted tmQCD continuum action, is given 

by 

Si= j da;"* {c iixj^^^F^,^^ + bf^fi^^pi; + ...}, (50) 

where the dots stand for further operators of dimension 5 (possibly includ- 
ing explicit mass factors) , which share the symmetries of the lattice action. 
The reason why I omitted them here is that they can be eliminated by the 
equations of motion. Furthermore, the second operator can be absorbed in 
an 0(a) shift of the standard quark mass parameter, so that one is really left 
with the SW term as the only relevant operator for on-shell improvement. 
Renormalised (connected) lattice correlation functions can be analysed in 
the effective theory, 

(O) = (O)™-^* - a(S'iO)™'^* + a((50)'=°"* + 0(a2), (51) 

where the cutoff dependence is explicit. We are here only interested in the 
leading cutoff effects at 0(a). To this order there are two contributions, 
first the insertion of the 0(a) part of the effective action ^i, and second 
the field specific counterterms SO. For example, with the choice 

= V;iix)P\y), (52) 

one finds the counterterm, 

SO = {cvid^Tl^^ix) + byti^Al{x)}P^{y) + ..., (53) 

where the dots stand for further terms which vanish by the equations of 
motion It should be emphasised that the 0(a) (and higher) corrections in 
the effective action are only treated as insertions, i.e. the expectation values 
^_^cont taken with respect to the continuum action 5*0. In writing down 
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the effective Symanzik theory there is thus an implicit assumption made, 
namely that one is working in the regime of continuum QCD where cutoff 
effects only appear as asymptotically small corrections. This assumption 
may certainly be wrong in some regions of parameter space, and particular 
care has to be taken in the presence of phase transitions. 
To proceed I introduce the 75T^-transformation, 

tjj ^ ijr^r'^ijjj ijj ^ ijji'-fc,T'^ , (54) 

which is part of the vector symmetry of two-flavour QCD. Hence 5*0 is 
invariant, but this is not the case for Si, i.e. one finds 

Sq — > Sq, Si — > — Si- (^^) 

For gauge invariant fields the transformation (54) squares to the identity, 
so that one may define an associated parity. For fields O with a definite 
75r^-parity one then finds, 

0^±0 ^ SO^tSO. (56) 

By applying the 75T^ transformation to the integration variables in the 
functional integral, one may derive identities between correlation functions, 
due to the invariance of the continuum action and functional measure. In 
particular, if we choose a 75r^-even field O, we find for the correlation 
functions at 0(a) 

(51 0)'^""* = -(51(9)^°"* =0, 
(50)™"* = -(50)™"* = 0, (57) 

and therefore 

(O) = (0)<=°"t+0(a2). (58) 
For a 75r^-odd O, one obtains 

which implies 

(O) = -a(S'iO)'^°"* + a{50y°''* + 0{a'^). (60) 

We may thus conclude that, at least in a small finite volume lattice corre- 
lation functions of 75T^-evcn fields are automatically 0{a) improved, while 
those of 75T^-odd fields vanish up to 0{a) terms. As a corollary, one may 



February 7, 2008 10:0 Proceedings Trim Size: 9in x 6in 



ms 



24 

state that standard Wilson quarks in a finite volume are automatically 
0(a) improved in the chiral limit. Although this is not the most interesting 
regime of QCD, it is somewhat surprising that this fact had not been no- 
ticed for more than 2 decades! To conclude this section note that in terms of 
the physical basis, (54) corresponds to the discrete flavour transformation, 

^' ^ -ir"^^', i>' i>'iT^. (61) 

A very similar argument based on parity transformations has been given by 
Shindler in ^. In a systematic analysis of the 75r^ symmetry (called Ti in 
this paper) can be found, showing that not only 0(a) but all odd powers of 
a vanish in 75T^-even correlators. This is not surprising, as this is implicit 
in the earlier analysis in ^, where the same conclusion was drawn. 

6.2.1. Uncertainty of the chiral limit 

If 0(a) improvement is automatic one might think that it should be possible 
to determine the critical mass m„ up to an intrinsic O(a^) uncertainty. This 
is not so, as I will now explain. The critical mass can be determined by 
tuning some PCAC mass to zero, and there is no obstacle for doing this 
in a finite volume. Now, the PCAC relation involves the axial current and 
density, A'^ and P", which have opposite 75r^-parities. According to the 
preceding discussion this means, for the first flavour components and with 
a 75T^-even source field Oeven, 

(5^4(a;)Oeven) = 2 mpcAC CP^)Oevc^ = 0(a2). (62) 
0(a) 0(a) 

The l.h.s. being 75T"'^-even must vanish up to O(a^), provided maximal 
twist is realised at least up to cutoff effects, i.e. rriR = 0(a). This implies 
that the PCAC mass is of 0(a), too, multiplying a correlation function 
which is 75T^-odd and therefore of 0(a). Thus no contradiction arises, the 
O(a^) of the l.h.s. is matched on the r.h.s. by two factors of 0(a). 

Another way to understand that an 0(a) shift in the critical mass does 
not ruin 0(a) improvement is to treat such a shift as an insertion of the 
standard mass operator ^■i/' into correlation functions. This operator is 
75T^-odd so that its insertion into a 75T^-even correlator produces an 0(o) 
effect, which together with the 0(a) mass shift yields an O(a^) effect. 
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6.3. Automatic 0(a) improvement in infinite volume 

When the infinite volume hmit is taken, the basic difference is the presence 
of spontaneous symmetry breaking and the appearance of non-analyticities 
in the mass parameters near the chiral hmit. As discussed ear her, twisted 
mass QCD is a valid regularisation of two-flavour QCD provided the con- 
tinuum limit is taken at fixed twist angle. To maintain maximal twist, 
i.e. a = 7r/2 one needs to tune the standard quark mass to mcr, which has 
an intrinsic 0(a) ambiguity. As long as the twisted mass is much larger 
than the typical 0(a) spread of Wcr, the twist angle may be considered 
well-defined, and the continuum limit is reached with O(a^) corrections. 
However, in practice one is interested in varying the quark mass at fixed 
cutoff, rather than studying the quark mass dependence only in the contin- 
uum limit. Approaching the chiral limit at fixed a by lowering the twisted 
mass one enters the regime where the twisted mass parameter becomes 
comparable to the 0(a) ambiguity of rricr- One may debate at this point 
whether the relevant comparison is with the uncertainty of mcr itself or 
rather with the size of typical O(a^) effects in correlation functions gener- 
ated by this uncertainty. In any case one reaches a point where the control 
over the twist angle is lost. When delivering my Nara lectures I interpreted 
this fact as a breakdown of the effective Symanzik theory. This is perhaps 
too rigid an interpretation. Rather one could say that for every definition of 
mcr, an effective twist angle is formed by the dynamics of the system, which 
may be far from the maximal twist one would like to maintain. Moreover, 
without further input is is impossible to know the effective twist angle for 
a given definition of Wcr- This is a disaster, as the whole interpretation 
of the theory rests on the twist angle, and a change in the effective twist 
angle (which remains unnoticed!) might strongly affect some correlators 
even at 0(1)! Fortunately this problem occurs close to the chiral limit, 
and thus in a region of parameter space where Chiral Perturbation Theory 
(xPT) is expected to describe the dynamics in terms of pion physics 
In particular, xPT is able to identify definitions of rricr in terms of pionic 
observables, which lead to an effective twist angle of a = tt/2, so that the 
Symanzik effective theory for maximally twisted mass QCD remains appli- 
cable in this region. For instance, this should be the case if one requires 
parity or flavour symmetry restauration, e.g. by imposing that a 75r^-odd 
pion correlation function vanishes. Note that the vanishing of the PCAC 
mass for a pion correlation function is a special case of such a condition. 
On the other hand, according to ^"^'^^ the condition of vanishing pion mass 
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(calculated in the untwisted theory) does indeed lead to a 0(1) variation of 
the effective twist angle. However, apart from larger cutoff effects of O(a^) 
this does not (yet?) seem to be a major problem in cf. figure 5. In 
any case, as the spontaneous symmetry breaking is closely related to the 
dynamics of pions, it seems that no statement can be made about generic 
definitions of rricr in a small volume, either from axial current conservation, 
or from parity or flavour symmetry restauration. 



0.55 p , I , I , I , I , I , I r 




0.07 



Figure 5. The continuum approach of Ftt in quenched tmQCD for various pion masses 
vs. a^/rp. 



7. Consequences of Parity and Flavour breaking 

The exact symmetries of lattice QCD with standard Wilson quarks include 
parity and flavour symmetry which are used to classify the hadron spec- 
trum. This is very convenient in any hadron analysis: even at fixed lattice 
spacing, the excited states which may occur in a given channel can be read 
from the Particle Data Book, with the exception of states with higher spin 
and/or angular momentum where the correspondence is spoilt by the lack 
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of rotational symmetry on the lattice. 

The situation is different in tmQCD since both parity and flavour sym- 
metry are broken by the Wilson term. As a consequence, the classification 
in isospin multiplets fails by terms of 0(a), or O(a^) if 0(a) improvement 
is at work. For instance, the neutral pion is not mass degenerate with the 
charged pions, or the nucleon A-resonances no longer form an exact isospin 
multiplet. Various simulations of quenched tmQCD have confirmed these 
expectations, and point to a restoration of flavour symmetry in the contin- 
uum limit 38,39,40^ although the expected rate oc a? for maximally twisted 
mass QCD has not in all cases been demonstrated convincingly. 

However, the splittings of isospin multiplets by cutoff effects are not the 
most serious drawback of parity and flavour symmetry breaking. In the 
spectral analysis of a hadronic two-point function all excited states with 
the same lattice quantum numbers may contribute. Even though the states 
violating continuum symmetries are multiplied by coefficients proportional 
to a, these states need to bo taken into account when working at fixed lat- 
tice spacing. Particularly annoying is the neutral pion, which shares all the 
lattice quantum numbers with the vacuum. One may thus add a neutral 
pion to any state without changing its lattice quantum numbers. The pres- 
ence of additional relatively light states may require a multistate analysis 
just to identify and subtract states which are a pure lattice artefact. More- 
over, correlation functions involving the light pion require the evaluation of 
disconnected diagrams. However, it should be emphasised that thesc^ prob- 
lems are purely technical; conceptually tmQCD is on a very solid basis, and 
in contrast to staggered fermions there is no mixing between flavour and 
spin degrees of freedom. 

7.1. N on- degenerate quarks and additional flavours 

Twisted mass QCD was originally formulated for a single doublet of mass 
degenerate flavours. This can easily be generalised to include more mass 
degenerate doublets. However, such a spectrum is quite unrealistic unless 

a non-degeneracy can be introduced within a doublet. Moreover, this non- 
degeneracy should not cause too much damage to all the nice properties of 
tmQCD. In particular, one needs to maintain the reality and positivity of 
the quark determinant, if such an action is to be used for simulations of 
full tmQCD. This is indeed possible, by introducing a mass splitting term 
as follows **^, 

C = -ii}{Ip+ m + in^^^T^ + 5mT^)ii), (63) 
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where Sm is the niass spUtting parameter. The mass spectrum is easily 
obtained by going to the physical basis and diagonahsing the mass matrix. 



Its eigenvalues are then found to be M± = ^m? + /Uq ± (5^. Translating 
this continuum situation to Wilson quarks in the obvious way, one first 
notices that the determinant of the twisted Wilson-Dirac operator must be 
real due to the conjugation property, 



Furthermore, the flavoTir structure of the determinant can again be reduced 
analytically, with the result, 



determinant at i5m = and continuity in 5m then imply positivity of this 
determinant for non-zero (5,,,. 

The mass splitting parameter is renormalised multiplicatively 5m,'K — 
Z^'^dm, where is the renormalisation constant of the non-singlet scalar 
density. As the positivity of the determinant follows from a condition on 
the bare parameters /iq and 5^, the corresponding condition in terms of 
the renormalised parameters involves a ratio of renormalisation constants, 
i-e. ^m,R < {Zp/Zs)iJ,R. The value of Zp/Zg depends on details of the reg- 
ularization, so that one cannot make a general statement about the ensuing 
limitations (if any). However, it is remarkable that one may use this ac- 
tion to perform numerical simulations with two non-degenerate light quark 
flavours, as needed for instance to study small isospin breaking effects. If 
used for strange and charm quarks 42,43^ however, one potentially has to 
deal with a fine tuning problem for the strange quark mass: for instance, 
assuming = 100 MeV and nic = 1300 MeV, these values are obtained 
as (700 ± 600) MeV. Finally, it should be said that the presence of the 
additional flavour non-diagonal breaking term renders the relationship to 
standard QCD more complicated, and the flavour structure needs to be 
dealt with explicitly in numerical calculations of quark propagators. 

8. A chiral twist to the QCD Schrodinger functional 

In order to solve scale dependent renormalisation problems the introduc- 
tion of an intermediate renormalisation scheme based on the Schrodinger 




(64) 
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functional (SF scheme) is an attractive possibility^^. Here I start by sum- 
marising its basic features in order to prepare the discussion of possible 
improvements . 

8.1. The QCD Schrddinger functional 

The QCD Schrodinger functional*^'^^(SF) is the functional integral for 

QCD where the Euclidean space-time manifold is taken to be a hyper cylin- 
der. The quantum fields are periodic in space, and Dirichlet conditions are 
imposed at (Euclidean) times xq = and xq = T. 

P+tp{x) Uo=o = P, P-tp{x) Uo=T= p', 
'4){x)P_ 1^0=0 = P, ■4>{x)P+ Uo=T= p', 

Ak{x) = Ck, Ak{x) U^=T= C'k, fc= 1,2,3, (66) 

with the projectors P± = ^(1 ±70). Correlation functions are then defined 
as usual, 

(0) = |z-i/" Oe-^j . (67) 

I ^fields J p=p'=0;p=p'=0 

O denotes some gauge invariant functional of the fields, possibly including 
the quark and antiquark boundary fields C, and C,, which are obtained by 
taking derivatives with respect to the quark boundary fields, viz. 

C(x).P-C(x) = ^, C(x).C(x)P- = -4^. (68) 

The name "Schrodinger functional" derives from the fact that such wave 
functionals arise naturally in the Schrodinger representation of Quantum 
Field Theory and the SF provides an example of a Quantum Field 
Theory defined on a manifold with a boundary. 

Using correlation functions derived from the Schrodinger functional, it 
is possible to define renormalised QCD parameters (the strong coupling and 
the quark masses), as well as renormalised composite operators (e.g. four- 
quark operators). Such renormalization schemes based on the Schrodinger 
functional (SF schemes) are attractive for the following reasons: 

• The finite volume is part of the scheme definition, i.e. all dimen- 
sionful quantities such as Euclidean time extent T, or boundary 

field parameters are scaled proportionally to L, the linear extent of 
the volume. As a consequence L remains the only scale in the sys- 
tem and can be identified with the renormalization scale by setting 
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H = L~^. Running parameters and operators then run with the size 
of the space-time volume, and one may apply recursive finite size 
techniques to bridge large scale differences (cf. subsect. 5.3.1) 

• SF schemes are made quark mass independent by imposing the 
renormalisation conditions in the chiral limit. Fortiuiatcly, the SF 
boundary conditions introduce a gap in the spectrum of the Dirac 
operator, which persists as the quark mass is taken to zero. This 
means that numerical simulations can be performed in the chiral 
limit, and no chiral extrapolation is needed to evaluate the renor- 
malisation conditions. 

• SF schemes are gauge invariant, no gauge fixing is needed. 

• Perturbation theory up to two loops is still feasable, due to the ex- 
istence of a unique absolute minimum of the action'^'*. This is to be 
contrasted with the situation on a hyper torus where perturbation 
theory becomes very intricate already at the one-loop level. 

• A further technical advantage consists in the possibility to use cor- 
relators involving zero momentum boundary quark and anti-quark 
fields. This is convenient in perturbation theory, and it leads to 
good numerical signals and reduced cutoff effects as compared to 
gauge invariant correlators in a periodic setting. 

All these nice properties come with a price: first of all, the presence of 
the boundary means that even the pure gauge theory suffers from 0(a) 
cutoff effects, caused by effective local operators of dimension 4, such as 
tr{FokFok} and tT:{FkiFki}, integrated over the boundary. When the quarks 
are included, there is even a dimension 3 operator, which can be absorbed 
in a multiplicative rescaling of the quark and antiquark boundary fields 
At order a, one expects dimension 4 operators like ip^oDoip and ip^kDk^ 
to contribute additional 0(a) effects ^. It is important to note that these 
cutoff effects are, unlike the 0(a) bulk effects of Wilson quarks, not due 
to the breaking of a continuum symmetry by the regularisation. Rather, 
such terms are to be expected with any regularisation of the Schrodinger 
functional. One may, however, write down a complete basis of 0(a) coun- 
tcrtcrms which contribute to a given observable. After reduction via the 
equations of motion, one typically ends up with 2-3 0(a) boundary coun- 
terterms. In practice it is then possible to monitor the size of the boundary 
0(a) effects by varying the coefficients. Perturbative results for these coef- 
ficients are often known to one-loop or even two- loop order and a non- 
perturbative determination may be conceivable. In summary, with some 
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extra work, the 0(a) boundary effects can be controlled and eventually 
eliminated. This is important, as otherwise the SF renormalisation proce- 
dure risks to introduce 0(a) effects even in 0(a) improved regularisations 
such as tmQCD at maximal twist or lattice QCD with Ginsparg- Wilson 
quarks. 



8.2. Decoupling of heavy quarks in SF schemes 

Quark mass independent schemes are very convenient to study the scale evo- 
lution for a theory with fixed quark flavour content. However, it also means 
that the decoupling of heavy quarks is not automatic, and one needs to 
match theories with different numbers of active flavours over quark thresh- 
olds. This is routinely done in perturbation theory, but it is not obvious 
that perturbation theory is adequate e.g. for matching the A^f = 4 and 
Nf = 3 effective theories over the charm quark threshold. One possibility 
to study decoupling consists in introducing a quark mass dependent SF 
scheme which would allow to study the non-perturbative evolution over the 
quark threshold until the heavy quark has decoupled. To define a mass 
dependent SF scheme it suffices to impose the renormalisation conditions 
at finite quark masses. Unfortunately, it turns out that the decoupling of 
a heavy quark in such a scheme is only linear in the inverse quark mass 
rather than quadratic. If the quark decouples very slowly, this means that 
it has to be kept longer in the evolution as an active degree of freedom, 
which could mean that widely different scales have to be accomodated on 
the same lattice. 

An example from perturbation theory^^'^'^ is given in figure 6. It shows 
the one-loop /9-function of the running coupling in the SF scheme as a 
function of z = mL, where m is some renormalised quark mass (its precise 
definition is not required to one- loop order). As z = mL is varied from to 
infinity, one expects to sec a smoothed out step function going from —1 to 
around the threshold z = 1. The solid and dotted curves (from 2 different 
SF schemes) do indeed show this behaviour, but the decoupling is rather 
slow compared to the MOM scheme^"'^ (dashed line). 

To understand this behaviour I propose a closer look at the Dirac oper- 
ator for free quarks and its spectrum in the continuum limit. 
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Figure 6. Decoupling of a heavy quark in the one-loop (3 function in two SF schemes 
and in the MOM scheme. See text for further explanation 



8.2.1. Free quarks with SF boundary conditions 

Let us consider a free quark tp in the continuum with homogeneous SF 
boundary conditions, 

P+V(a;) U„=o= 0, P-i>{x) U„=T= 0. (69) 

Then 75 (^+m) is a hermitian operator with smooth eigenfunctions and no 
zero modes Evaluating the eigenvalue equation for any of its eigenfunc- 
tions (p at the boundaries one finds, 

P+75(^+to)(/? |a;o=o= ^ {do - m)P-'f \x„=o= 0, 
P-j5{P+m)>fiU,=T=0 => (5o + m)P+^ U=T= 0. (70) 

The complementary components thus satisfy Neumann conditions mod- 
ified by the mass term m. The eigenvalues A are of the form A = 
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±-\/pq + + TV?, where po is determined as non- vanishing solution of 

tan(poT') = ^Po/™- It is obvious from this equation that ^-nd thus A 
are not symmetric under m ^ —m. This is generic and can be understood 
as a consequence of chiral symmetry breaking by the boundary conditions. 
As a result one expects, for any observable in the SF, the asymptotic small 
mass behaviour cx m (rather than m,^), and similarly for heavy quarks the 
corrections oc 1/m (instead of 1/m^), as illustrated in figure 6. 

At least for even numbers of flavours a possible way out consists in 
adding a twisted mass term and setting m = 0. Then j^t^(P+ ifiq'y^T^) 
is again hermitian. With this Dirac operator, the complementary field 
components at the boundaries satisfy simple Neumann conditions, and the 
spectrum is symmetric under a change of sign of the twisted mass term. 
A physically equivalent solution is obtained by staying with the standard 
mass term and rotating the boundary projectors instead. This will be 
discussed in more detail below. However, a caveat remains as only the 
simultaneous decoupling of an even number of quarks can be studied in 
this formalism. On the other hand, it may be sufficient to compare to 
perturbative decoupling in this slightly unphysical setting, in particular if 
a perturbative treatment turns out to be satisfactory. 

8.3. SF boundary conditions and chiral rotations 

Let us consider flavour doublets tp' and tjy which satisfy homogeneous stan- 
dard SF boundary conditions. Performing a chiral rotation, 

■0' = cxp(m75T'^/2)'0, = -0 exp(m75r'''/2), (71) 

one finds that the fields tp and 'ip satisfy the chirally rotated boundary 
conditions, 

P+(a)0(x) U„=o= 0, P_(a)^(,T) U.„^T= 0, 
i;{xhoP^{a)\^,=o=0, i;{x)-foP+{a) \^,=t= 0, (72) 
with the projectors, 

P±(a) = i[l ±7oexp(^a75r3)]. (73) 
Special cases arc a = and a = Tr/2 where one obtains, 

P±{0) - P±. ^±(V2) = 0± = i(l ± no75T3). (74) 

We perform again a change of variables in the functional integral. In- 
cluding mass terms as well, we label correlation functions by a subscript 
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{m, iiq, P+{a)), i.e. we include the projector defining the Dirichlet compo- 
nent of the quark field at xq = 0. The generalisation of formula (22) then 
reads: 

{0[tp,i}]){m,,,^,p±) = (0[i?(a)V', V'^(Q;)]>(TO,/i„p±(a)), (75) 
with mass parameters m and flq given by 

rh = m cos a — fiq sin a, flq = m sin a + fiq cos a. (76) 

The boundary quark fields arc included in this transformation by replacing 

C(x) ^ V5(0,x)P+, C(x) ^ P-^(0,x). (77) 

This extends the equivalence between correlation functions of tmQCD and 
standard QCD to correlation functions derived from the Schrodinger func- 
tional. Simple examples are provided by purely gluonic observables 0[U], 
such as the SF coupling constant. Eq. (75) then implies, 

= {om (78) 

In other words, either the mass term is twisted and one stays with standard 
SF boundary conditions, or the mass term is standard and the boundary 
conditions are fully twisted. In both cases one expects a quadratic de- 
pendence on the mass parameter and hence a relatively fast decoupling of 
heavy quarks. 

8.4. SF schemes with Wilson quarks and 0{a) improvement 

From the discussion of 0(a) improvement in section 6 one may conclude 
that 75T^-even observables computed with Wilson quarks in a finite volume 
and with periodic boundary conditions are automatically 0(a) improved at 
zero quark mass. As SF schemes are usually defined at zero quark mass, 
it seems natural to ask how the SF boundary conditions interfere with this 
property. It is useful to think of 0(a) effects to arise from different sources. 
First there are the 0(a) boundary effects, which are cancelled by intro- 
ducing the 0(a) boundary counterterms to the action and the boundary 
quark and antiquark fields. Second there are 0(a) effects from the bulk 
action which may be cancelled by the Sheikholcslami-Wohlcrt term, and 
third there are the 0(a) effects associated with the composite operators 
in a given correlation function. It is interesting to note that 0(a) cutoff 
effects from the bulk action are often quite large in SF correlation func- 
tions. This is illustrated in figure 7 which shows the relative cutoff effects 
in the perturbative one-loop coefficient of the step-scaling function of the 



February 7, 2008 10:0 Proceedings Trim Size: 9in x 6in 



ms 



35 

four-quark operator needed for Bk The operator here is unimproved, 
and the boundary effects remain uncancelled in order to mimick the non- 
perturbative procedure of The most dramatic reduction of cutoff effects 
occurs when the Sheikholeslami-Wohlert term is included. Moreover, this 
has the side effect to reduce the ambiguity in the zero mass point, so that 
with the standard SF it makes sense to implement 0(a) improvement even 
if it is not complete. 



• SW, m_^{L/a) _ 
■ Wilson, mjL/e.) - 
i SW, m^(=») 

• Wilson, m^(">) 




a/L 



Figure 7. Relative cutoff effects in the one-loop coefficient of the step-scaling function 
of the Bji operator. Shown are two different regularisations (csw = 0, 1) with two 
definitions of the zero mass point. 



8.5. The Schrodinger functional and 0{a) improvement 

The reason why automatic 0{a) improvement fails is that the 75T^- 
transformation (54) changes the projectors of the quark boundary con- 
ditions, 

P±15t' =-f5r'P^. (79) 
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The boundary conditions, just like mass terms, define a direction in chiral 
flavour space. Tliis means tfiat tfie 75T^-transformation yields inequivalent 
correlation functions even in the chiral limit. For a 75 -even operator O 
one finds 

{0){m,ti^,P+) {0)(-m,n^,P-)- (80) 

It thus appears that the standard SF does not allow for the definition of 
75T^-even correlation functions, and bulk 0(a) improvement is not auto- 
matic. 

A possible solution is obtained by changing the projectors used to specify 
the Dirichlet components such that they commute with j^t^. Allowing for 
an additional flavour structure one may think of ^(1 ± 7ot^) or 

g± = i(l±i7o75T-''). (81) 

Interestingly, the projectors Q± also appear in the chiral rotation of the 
SF by a = 7r/2. Besides automatic 0(a) improvement, the implementation 
of these boundary conditions may lead to some interesting checks of uni- 
versality by comparing SF correlation functions in the standard framework 
and at maximal twist. Note that this direct comparison was not possible in 
where a twisted mass term was introduced whilst keeping the standard 
SF boundary conditions. 

8.6. The SF with chirally rotated boundary conditions 

The implementation of some given boundary conditions is not straightfor- 
ward on the lattice, and some care has to be taken to ensure that one really 
ends up with the desired continuum theory. A successful implementation 
of the maximally twisted boundary conditions involving the projectors Q± 
has been described in and relies on an orbifold construction to ensure 
the correct continuum limit. 

8.6.1. Symmetries and Counterterms 

Apart from the absence of a dimensionful parameter, the symmetries of the 
SF with maximally twisted boundary conditions are identical to those of 
tmQCD. One may then list the possible boundary counterterms of dimen- 
sion 3 allowed by the symmetries: 

Ki = 4>i^^T^i^, K± = ^Q±^p. (82) 

As time reflection combined with a flavour permutation is a symmetry of 
the SF, it is enough to discuss the counterterms at a;o = 0. Ki corresponds 
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to the logarithmically divergent boundary counterterm in the standard SF, 
which leads to a multiplicative renormalization of the quark boundary fields. 
The operator Kj^ only involves Dirichlet components at xq = and is 
therefore irrelevant for most correlation functions used in practice. The 
remaining operator K_ only contains non-Dirichlct boundary components. 
If rotated back to the primed basis it becomes proportional to ^'i^5T''^P-tp' , 
which violates flavour symmetry and parity just like a twisted mass term. 
As these are symmetries which are restored in the continuum limit one 
concludes that this counterterm must be scale-independent. Its coefficient 
can be fixed by requiring that a parity violating SF correlation function 
vanishes at finite a. 

This analysis can be extended to dimension 4 operators which appear 
as 0(a) boundary counterterms It turns out that the situation is com- 
parable to the standard SF, i.e. there are a couple of counterterms which 
one needs to tune in order to eliminate the 0(a) boundary artefacts. 

8.7. An example from perturbation theory 

In perturbation theory, the values of all boundary counterterms are known, 
so that one may study both the equality of properly matched standard 
and twisted SF correlation functions, and confirm automatic bulk 0(a) 
improvement. A first example is given by the SF coupling, which can be 
related perturbativcly to the MS-coupling, 

f(L) - g^ifi) + k,UiL)g^{fi) + 0{g'). (83) 

The fermionic contribution to the one-loop coefficient, fci = fci,o + A^ffci,!, 
has been computed in yielding fci^i = — 0.039863(2)/(47r). In practice, 
the perturbative data is obtained for a sequence of lattices, and one then 
expects the asymptotic large L/a behaviour: 

f{L/a) ~ ro + (a/L)[ri + si ln(a/L)] + 0{a^). (84) 

Here vq = ki i is the continuum limit value, and the 0(a) effects lead to 
non- vanishing values of ri and si. In the standard Schrodinger functional 
set-up one expects that ri is eliminated by the boundary counterterm pro- 
portional to tr(Fofc-F'ofc), whereas si is due to bulk 0(a) effects from the 
action, and thus proportional to ciw — 1- On the other hand, with twisted 
SF boundary conditions one expects that ro remains the same, due to uni- 
versality, ri is cancelled again by a boundary counterterm, and si should 
vanish independently of the value of Csw-This expectation is indeed con- 
firmed numerically. A similar test can be performed with the tree level 
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quark propagator in a non-vanishing gauge background field, induced by 
choosing non- vanishing gauge field boundary values and C^. One then 
expects that, with the correct tree-level boundary counterterms, the bulk 
0(a) lattice artefacts will again be either proportional to Csw — 1 (standard 
SF) or absent (twisted SF). Again this expectation is confirmed. However, 
in contrast to the SF coupling this test can not be extended beyond the 
tree level, unless one fixes the gauge. 



9. Conclusions 

Lattice QCD with Wilson type quarks remains an attractive regularisation 
of lattice QCD. Some of its problems can be alleviated by introducing a 
chirally twisted quark mass term. While the theories remain equivalent in 
the continuum limit, the twisted mass term supplies an infrared bound on 
the spectrum of the Wilson-Dirac operator which renders the quenched and 
partially quenched approximations well-defined. Some of the notorious lat- 
tice renormalisation problems of standard Wilson quarks can be by-passed, 
and tmQCD at maximal twist is automatically 0(a) improved. These ad- 
vantages are balanced by parity and flavour breaking and the fact that 
tmQCD comes naturally with an even number of quarks. 

The Schrodinger functional has become an indispensable tool to tackle 
non-perturbative renormalisation problems in lattice QCD However, the 
standard set-up leads to a slow decoupling of heavy quarks, and is in con- 
flict with automatic 0(a) improvement of masslcss Wilson quarks. This 
motivates the application of a chiral twist to the SF boundary conditions. 
It is thus possible to extend equivalence between tmQCD and standard 
QCD to correlation functions derived from the Schrodinger functional. This 
allows for interesting tests of universality and the maximally twisted SF is 
compatible with automatic 0(a) improvement, as I have illustrated with 
simple perturbative examples. 
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